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Several recent studies have been performed in the Icing Research Tunnel (IRT) at NASA 
Glenn Research Center focusing on the evolution, spatial variations, and proper scaling of 
ice roughness on airfoils without sweep exposed to icing conditions employed in classical 
roughness studies. For this study, experiments were performed in the IRT to investigate the 
ice roughness and thickness evolution on a 91.44-cm (36-in.) chord NACA 0012 airfoil, swept 
at 30-deg with 0° angle of attack, and exposed to both Appendix C and Appendix O (SLD) 
icing conditions. The ice accretion event times used in the study were less than the time 
required to form substantially three-dimensional structures, such as scallops, on the airfoil 
surface. Following each ice accretion event, the iced airfoils were scanned using a ROMER 
Absolute Arm laser-scanning system. The resulting point clouds were then analyzed using 
the self-organizing map approach of McClain and Kreeger to determine the spatial 
roughness variations along the surfaces of the iced airfoils. The resulting measurements 
demonstrate linearly increasing roughness and thickness parameters with ice accretion time. 
Further, when compared to dimensionless or scaled results from unswept airfoil 
investigations, the results of this investigation indicate that the mechanisms for early stage 
roughness and thickness formation on swept wings are similar to those for unswept wings. 


Nomenclature 


= accumulation parameter 
angle of attack 
airfoil maximum roughness, the maximum of the RMH values along an airfoil or wing surface 
codebook vectors 
= Acoefficient related to the transient scaling function (~ 0.5 for unswept wings) 
= neighborhood function of i to j codebook vectors 
= the number of surface points in the neighborhood of a specific SOM codebook vector 
= codebook vector index 
= droplet inertial parameter 
= Length scale related to the Stokes number 
= liquid water content [gm/m*] 
= median volumetric diameter [um] 
= airfoil or mean ice shape surface normal coordinate direction 


the orthogonal distance from a surface point (x/) to the mean surface manifold as described by the 
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SOM codebook vectors 


Nor = Fully-dense Rime ice stagnation thickness on unswept wings 

Nar = Fully-dense Rime ice surface-orthogonal thickness on swept wings 

N, 0 = the measured ice orthogonal-thickness at each codebook vector relative to the clean wing surface 
N, 0.AL = the measured attachment line orthogonal-thickness relative to the clean airfoil or wing surface 
n = stagnation point freezing fraction 

RMH = 99%-Gaussian roughness maximum height evaluated at each codebook vector (=3.09R,) 

Rg = the root-mean-square or “standard deviation” roughness height 

Ya = leading edge radius of curvature 

SOM = Self-Organizing Map 

S = mean ice shape surface tangential coordinate direction 

So = the clean airfoil tangential coordinate direction and distance 

S(c) = Spatial scaling function 

Sk = Stokes number 

3 (z) = Transient (time-based) scaling function 

T static = the freestream static temperature 

Tiotal = the freestream stagnation or total temperature 

VorUx = Freestream velocity 

x = element of point cloud data set 

x = arandom variable 

a = local direction angle of manifold through a codebook vector or airfoil angle of attack 

hp = local surface collection efficiency at the location of the liquid-film breakdown 

Bs = local surface collection efficiency 

Bo = straight wing stagnation collection efficiency 

Ba = swept wing attachment-line collection efficiency 

y = direction angle of surface point relative to manifold direction through winning codebook vector 
Vs = direction angle of a surface point relative to the airfoil design coordinates 

Yi = direction angle of the surface point of liquid-film break down relative to the airfoil design coordinates 
At; = the ice accretion event time 

A = Sweep angle 

Pice = density of ice 

Pw = density of liquid water 

o = Scaled position parameter 

T = Scaled time parameter 

Hair = viscosity of air 


I. Introduction 


OUGHNESS on iced airfoils increases the skin friction drag, increases the rates of heat transfer from the 

freezing water layer to the airstream, and induces higher levels of boundary-layer turbulence than would be 
observed on mean representative ice shapes without roughness. Further, ice roughness increases the droplet capture 
efficiency of the local surface relative to a smooth surface. Because it couples the fluid flow, heat transfer, and 
droplet impingement processes, ice roughness formed during the early stages of the ice accretion process is an 
important factor in the overall in-flight ice accretion process and is thought to ultimately affect the resulting ice 
shapes exhibited much later in the ice accretion process. 

Because of its importance to the overall ice accretion process, the study of roughness on ice accretion surfaces 
has a long history. However, most of the historical studies have involved image analysis approaches to characterize 
roughness element morphology. References [1-3] represent examples of the morphological roughness investigations. 
Following the development of laser scanning techniques in the Icing Research Tunnel (IRT) at NASA Glenn 
Research Center [4-7], several investigations of ice roughness evolution [8,9] and ice roughness scaling [10] have 
been performed on straight wings with symmetric (NACA 0012) airfoil shapes. More specifically, the roughness 
evolution investigations of McClain et al. [8] and McClain et al. [9] focused on the temporal evolution of roughness, 
while the scaling investigation of McClain et al. [8] focused on the spatial variations of roughness features by 
identifying the proper scaling quantities relative to the stagnation point cloud collection efficiency. 
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While the previous laser-scanning studies have provided significant insights into roughness temporal and spatial 
evolution on straight wings, roughness evolution on swept wings have yet to be addressed. Vargas [11] presented 
review of swept-wing roughness mechanisms. Vargas [11] indicated a need to better understand swept-wing 
roughness evolution and morphology because of their importance on scallop and lobster-tail formation on swept 
wings. Because the boundary-layers on swept wings are three-dimensional, the morphology of the roughness 
elements are thought to be critical to promoting instabilities in the boundary-layers and leading to the formation of 
the highly three-dimensional ice features formed and observed on swept wings. 

The objectives of this study were to measure ice roughness and thickness topographies on swept wings during 
short-time icing events, that is, before the formation of large three-dimensional features such as scallops or lobster- 
tails. For the study, the same basic Appendix C [12] and SLD [13] conditions of McClain et al. [9] and McClain et 
al. [8], which were based on the classical investigation of Anderson et al. [3], were employed. The temporal 
evolution and the scaled spatial evolution results will be compared to the past studies on straight wings to investigate 
similarities between straight and swept-wing roughness evolution. By illuminating the physics of swept-wing 
roughness evolution, the effort seeks to improve the heat transfer model in NASA ice accretion codes, LEWICE 2D 
and LEWICE 3D, which will be used in the design of current and future air vehicles employing swept wings or 
blended surfaces. 


II. Roughness Scaling Concepts 


The prior straight-wing roughness and evolution studies [8-10] have approached the issue of roughness temporal 
evolution and spatial variations as separable aspects of the roughness formation process. That is, the temporal 
evolution and the spatial variations along the airfoil surface are modeled using the product of a scaled time function 
and spatial function. 


RMH (t, Sy, LWC, MVD...)= 3(z)5(c) (1) 


In Eq. (1), 3 (c) is a function that depends on the icing event time and the supercooled cloud properties related to the 
impingement process, and S(c) is function that relates to the surface position, that is its geometry relative to the 


upstream flow, and to the supercooled cloud properties related to the liquid film dynamics. 
The temporal function relates to the exposure time of the aircraft surface in a supercooled cloud. For traditional 
aerosol-and-air problems, the Stokes number is defined as 


_ PyU,MVD* 
18h gi-L 


Sk (2) 


Following the development of Tsao and Lee [14] and Tsao and Kreeger [15], when the length scale employed in the 
Stokes number definition is the leading edge radius of the airfoil, the result is commonly referred to in the icing 
literature as the droplet inertial parameter, K: 


= PU ..MVD* 
1B Maia 


K (3) 


To maintain droplet trajectory similarity, the Langmuir and Blodgett [16] stagnation point collection efficiency, Bo, 


is employed where 
1 
Loe") 


114 Ko -3) 


In Eq. (4), Ko is the modified inertial parameter of Langmuir and Blodgett [16], defined as 


Bo = (4) 
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eee 5 


where the parameter A/A,, is the droplet range parameter and is defined as a function of the droplet Reynolds 
number, Res = pV MVD/ 1, as 


A 1 
7 (6) 
Asx 0.8388 + 0.001483 Re; + 0.01847,/Re 5 


Once the stagnation point collection efficiency is known, the volume of ice impinging an incremental area, as 
depicted in Figure 1, may be evaluated as 


_ LWC-V- By - At, 
Pice 


d¥, 


ice 


dA (7) 
Dividing by the incremental area and multiplying by 270/2ro results in what McClain et al. [10] referred to as the 


leading-edge, fully-dense, theoretical rime ice (n = 1) thickness, No,r, for straight wings and airfoils. 


_ LWC-V- By At, 219 
Pice 2K 


0,R = 2mA, fo (8) 


In Eq. (8), Ac is the accumulation parameter defined as 


LWC-V-At 
7 eens 


9 

2K 0° Pice ad 
For straight wings, the prior laser-scan based scaling investigation [10] expressed the temporal function in terms of 
the 99% maximum roughness height was found to be approximately 


3(r) = AMR(At, , cloud properties ) =C,Nor (10) 


Revisiting McClain et al. [10], the constant C; for unswept wings was found to be approximately 0.5, and captured 
both the roughness evolution versus time and the stagnation point ice thickness versus time. 

For swept wings with surfaces oblique to the flow, dA in Figure 1 is the projected area in the direction of the 
flow, and the Rime ice would be expected to spread over the actual surface area. Consequently, the fully-dense, 
theoretical Rime ice thickness definition of Eq. (8) would be modified as shown in Eq. (11) to represent the swept 
wing leading-edge “surface-orthogonal” Rime ice thickness, Na,r. 


Nar = 2%A-By cos(A) (11) 
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Figure 1. Supercooled Cloud Impinging an Incremental Area on a Stagnation Line 
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For straight wings, McClain et al. [10] found that surface roughness variations along the ice surface were 
separated into two regions, as demonstrated in Figure 2(a), and defined based on the local collection efficiency 
corrected by the local surface angle relative to the flow. In the first region, the roughness features were suspected as 
being formed by the liquid film instability dynamics, while in the second region, collection efficiency dynamics 
governed the roughness variations. While an expression was not provided by McClain et al. [10], near and aft of the 
glaze ice plateau where the surface collection physics are dominant, the surface variations are reasonably well 
described using Eq. (12), which is represented by the gray, dashed curve in Figure 2(a). 


5(o)=__" * = 4 
C.Nor Weil0.875(B p, — 0.15} B y, —0.15,3) 


ne RMH(8,. 2 5) WeilB, —0.05, B , —0.15,3) @ 


as tanh[10(6 », - 8, ) (12) 


Where 
Bs = Bo cos(y, +a) (13) 
Bp = Bo cosly p +a) (14) 


and 7 is the angle of the surface relative to the airfoil design flow coordinates, as depicted in Figure 3, @ is the angle 
of attack, and 7 is surface angle relative to the airfoil design coordinates at the location of the liquid-film break 
down and ultimately the location of the glaze ice plateau. Further, Wei(x,/,k) in Eq. (12) is the Weibull probability 
density function described as 


" k-l : 
Wei(x,A,k) = *(=) e iA) (15) 


Where x is a positive real number, A is referred to as the scale parameter, and k is called the shape parameter or often 
referred to as the Weibull modulus. 

In equation (12), C; is set to 0.5, and the local collection efficiency at the point of film breakdown (£#) was set to 
0.45. These values are represented in Figure 2 by the black dashed lines. In both Figure 2(a), which presents the 
scaled roughness heights for the unswept wing results, and Figure 2(b), which presents the scaled ice thickness 
measurements for the unswept wing results of McClain et al. [10], the horizontal dashed lines represent a C; value of 
0.5 based on the ordinate variable definition and an ordinate value of 1.0. The collection efficiency at the point of 
film breakdown (fy) is identified in both Figure 2(a) and 2(b) by the vertical dashed line in each subfigure. In 
Figure 2(b), the location of fw is more easily identified as the location where the thickness measurements transition 
from the stagnation region to the collection efficiency region, which was also identified by McClain et al. [10] as the 
edge of the glaze-ice plateau. 

For swept wings, roughness in regions where the collection efficiency physics are dominant would be expected 
to follow a variation similar to Eq. (12). However, the local collection efficiencies are expected to require alteration 
as shown in Eqns. (16) and (17) based on the wing sweep angle. 


B,=Bo cos(y, +a)cos(A) (16) 
Bw» = Po cos|y g, +a)cos(A) (17) 


The experiments performed were constructed to determine if the modifications to the scaling approach for wing 
sweep in Eqns. (11), (16), and (17) are sufficient to capture surface roughness variations and temporal evolution on 
swept wings. If the temporal and spatial modifications capture the development of roughness on swept wings, then 
the results would indicate that the underlying physical mechanisms for roughness evolution on swept wings are the 
same as the roughness formation mechanisms on straight wings. 
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Figure 2. Scaled Roughness Height and Thickness Variations versus Local Collection Efficiency for Unswept 


Wings [10]: (a) Roughness and (b) Thickness 
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Figure 3. Unswept Airfoil Geometry Relating to the Local Surface Collection Efficiency (f) Relative to the 
Stagnation Point Collection Efficiency (fo) 


Ill. Methodology 


To investigate ice roughness and thickness evolution on a swept airfoil and their similarities to roughness and 
thickness evolution on straight airfoils, the ice shapes were 1) created, 2) measured, and then 3) characterized using 
the SOM approach of McClain and Kreeger [17] for a wing swept at 30°. The following sections describe each of 
the ice roughness generation and characterization steps. 


A. Ice Shape Generation 

All of the experimental measurements were performed in the IRT over two nights of testing. A 91.44-cm (36- 
in.) chord NACA 0012 airfoil, swept at 30°, was used for the study. The swept airfoil spanned 152.4-cm (60-in) 
long. One of the nights of testing was devoted to Appendix C conditions, and the other night was devoted to 
Appendix O or supercooled large droplet (SLD) icing conditions. 

The flow conditions chosen for the Appendix C tests were essentially the same conditions employed by 
Anderson et al. [3] in the historical studies of ice roughness on unswept wings with a fixed MVD of 30um. Table 1 
presents the freestream static and total temperatures, the freestream velocity, the median volumetric diameter, the 
liquid water content, the accumulation parameter, the freezing fraction, and the stagnation point collection 
efficiency. As shown in Table 1, the primary experimental dimensions varied during the tests as the ice accretion 
time resulting in accumulation parameters ranging from 0.03 to just under 0.5. 

Table 2 presents the cases for the SLD tests. All of the SLD tests were performed using an MVD of 150 um. 
Comparing Tables | and 2, the static air temperatures and LWCs were kept consistent between all tests. However, 
the freestream velocity for the SLD cases was slightly lower than for the Appendix C cases. Combined with the 
larger droplet size, the slightly lower velocity, which is still within the IRT velocity controllability of +0.5 m/s [18], 
resulted in lower stagnation point freezing fractions, n, for the SLD cases (0.19 versus 0.22 for the Appendix C 
cases). 

For each test, the airspeed and freestream total temperature were set and the spray bar air and water pressures 
were selected to provide the appropriate LWC, MVD, and freezing fraction. A thermocouple embedded in the 
NACA 0012 airfoil was used to determine when the airfoil had reached thermal equilibrium with the flow. Once the 
thermocouple reported the static temperature of the flow indicating that it was in thermal equilibrium with the flow, 
the spray bars were actuated and closed after the predetermined spray time. 


B. Ice Shape Measurement 

Following the completion of the icing spray, the wind tunnel velocity was reduced to less than 10 knots (5 m/s) 
while keeping the static temperature around -4 °C to avoid thawing of the ice shape. The iced airfoil was painted 
using a tetrahydrofuran (THF)-based titanium dioxide paint. A ROMER Absolute Arm laser scanning system was 
placed in the IRT test section upstream of the airfoil. A scan was then made of the leading 200-250 mm of the 
airfoil leading edge (in the streamwise direction) on both sides of the airfoil downstream of the leading edge. The 
scans were approximately 150-mm wide (in the spanwise direction) and were performed at the location of the airfoil 
corresponding to the center of the test section. 
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Table 1. Summary of Test Parameters for the Appendix C Tests 
Case T static Ttotat MVD LWC Ats 


Number | Cc) | CC) | YS) | Gum | cm) | Gee | * 7 fo Ba 
012816.01 | 46 | -23 67.1 30.0 0.60 20 | 0.030 | 022 | 0.713 | 0617 
012816.02 |_-4.6 | -23 al 30.0 0.60 40 |_0.060 | 0.22 | 0.713 | 0.617 
012816.03 | 4.6 | -23 67.1 30.0 0.60 60_|_0.091 | 022 | 0.713 | 0617 
012816.04 | 4.6 | -23 67.1 30.0 0.60 80/0121 | 0.22 | 0.713 | 0617 


012816.05 | -4.6 -2.3 67.1 30.0 0.60 100 0.151 0.22 0.713 0.617 
012816.06 | -4.6 -2.3 67.1 30.0 0.60 160 0.242 0.22 0.713 0.617 
012816.07 | -4.6 -2.3 67.1 30.0 0.60 240 0.363 0.22 0.713 0.617 
012816.08 | -4.6 -2.3 67.1 30.0 0.60 320 0.484 0.22 0.713 0.617 


Table 2. Summary of Test Parameters for the SLD Tests 


Case T static Ttotal MVD LWC Mts 

Number | eC) | cco) | Y™) | Gm | (emmy | (ec) | 4° ie fo Ba 

01291601 | 4.6 | 24 | 66.7 | 150.0 | 0.60 30 10.030 | 0.19 | 0.953) 0.825 
012916.02 | 4.6 124 166.71 150.0 10.60 | 40 10.060 [0.19 | 0.953 | 0.825 
012916.03 | 4.6 124 | 66.7 1150.0 | 0.60 60 10.090 10.19 | 0.953 | 0.825 
012916.04 | 46 |..24 166.7} 150.0 | 0.60 80 10.120 10.19 | 0.953 | 0.825 
012916.05 | 4.6 124 166.71 150.0 10.60 | 100 10.150 [0.19 | 0.953 | 0.825 
012916.06 | 4.6 124 166.71 1500 | 0.60 | 16010241 [0.19 [0.953 | 0.825 
012916.07 | 4.6124 166.71 150.0 | 0.60 | 24010361 [0.19 | 0.953 | 0.825 
012916.08 | 4.6 124 | 66.71 150.0 10.60 | 320 10481 [0.19 [0.953 | 0.825 


C. SOM Roughness and Thickness Evaluation 

The surface point clouds were then processed using the Self-Organizing Map (SOM) approach developed by 
McClain and Kreeger [17]. The self-organizing map, or sometimes referred to as a Kohonen Map, is a clustering 
method for the detection of non-linear manifolds, which may be curves or surfaces, in multi-dimensional space 
[19]. SOMs depend on the use of codebook vectors, b, which may also be called codebook points or neurons, to 
represent clumps of data. Following convergence of the SOM method, each codebook vector will be located at the 
spatial centroid of the clump of data that it represents. In its simplest essence, self-organizing maps are employed to 
capture trends of large data sets by representing those large data sets by a relatively small set of codebook vectors. 
When applied to an ice shape, the intent of performing the SOM is to extract the “form” of the ice shape from the 
surface variations. For detailed information on self-organizing maps and their application for iced airfoil description 
and roughness evaluation, please consult Refs. [17] and [20]. 

When applied to an iced airfoil point cloud without sweep or significant spanwise shape changes, the SOM is 
expected to identify a curve in the Chord-Chord Normal plane (or x-y plane when the spanwise axis is placed in the 
z-direction), which represents the mean shape of the rough airfoil. The nature of the SOM method and the 
positioning of the codebook vectors along a “daisy-chain” enable a statistical evaluation of iced airfoil surface 
roughness. Since the “clumps” of points are distributed about the codebook vectors, the deviations of the point 
measurements in the clumps can be used to evaluate the coverage statistics and uncertainty of the codebook vector 
representation. 

Figure 4 shows a single surface measurement, x/, and its closest codebook vector b". The two neighboring 
codebook vectors along the daisy-chain of codebook vectors representing the manifold are also shown. In the 
approach used for this study, the manifold is assumed to be a first-order manifold in two-dimensional space with the 
characteristic that at each codebook vector the local slope of the manifold is equal to the central finite-difference 
evaluated using the two closest surrounding codebook vectors. The approach used assumes that all deviations from 
the manifold are normal to the manifold. That is, the deviation of a surface measurement normal to the line through 
the codebook vector with the local slope set by the neighboring codebook vectors is considered the “height” of the 
surface point above or below the local manifold. 
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Figure 4. Metrics of Local Point about a Codebook Vector [9] 


Based on the SOM-manifold description, the root-mean-square roughness height is calculated at each codebook 
vector as 


iw ; 2 
R,= SNe (17) 
j=l 


In some roughness studies, the 99.9% roughness maximum height (RMH) based on a Gaussian distribution is 
calculated using 3.09 times the root-mean-square roughness height. 


1 
iw ; 2 
RMH ,, =3.09R, =3.09 ae (18) 
j=l 


Where J is the number of surface points for which b” is the winning (closest) codebook vector. The RMH is the 
local 99.9%-maximum roughness height and is evaluated at each codebook vector based on the local or 
neighborhood statistics. The RMH will vary along the surface arc length. While the RMH is not a universal 
roughness descriptor, the RMH value is important for this study. Since the RMH represents a Gaussian prediction 
of the 99.9% maximum distance from the mean elevation to the tallest peaks in a data set, the RMH is a reasonable 
tool to compare the statistical results to the morphological descriptors such as roughness element diameter and 
height used in the historical roughness studies. More details regarding the roughness evaluation and the associated 
measurement uncertainties may be found in McClain [19]. 

The ice thickness relative to the clean airfoil shape, No, is calculated in a fashion very similar to the ice 
roughness height, as demonstrated in Figure 4, with the exception of instead of using the codebook vector as the 
surface reference, the original airfoil design coordinates are used. Further details regarding the ice thickness 
evaluation may be found in McClain [19]. 

The point clouds were analyzed using the Airfoil Roughness Evaluation System (ARES), which was developed 
to automate the roughness characterization approach of McClain and Kreeger [17]. The approach of McClain and 
Kreeger [17] was originally implemented using multiple software tools. ARES is a set of Matlab functions that 
performs the SOM characterization of the mean ice shape and then performs the two-dimensional statistical analyses 
required to characterize the roughness variations along the airfoil surface in the flow direction. After being read into 
ARES, the point clouds were rotated to be aligned with the primary axis of the wing. Consequently, the roughness 
and thickness measurements are performed in the projected X-Y airfoil coordinate system (as opposed to the wind 
tunnel coordinate system) and result in “surface orthogonal” roughness and thickness measurements. 
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IV. Results and Discussion 


The results of the investigation are divided into three separate aspects: 1) visualization of roughness and 
thickness topographical variations, 2) temporal roughness and thickness evolution, and 3) spatial roughness 
variations and scaling. Each aspect of the results is described in the following subsections. 


A. Visualization of Roughness and Thickness Topography 

Figure 5 presents a visual comparison of the roughness and thickness topographies for the 012816.05 case 
compared to the 012916.05 case. Both cases have the same icing event times, At, , of 100 seconds. However, the 
012816.05 case is an Appendix C case (MVD = 30m), and the 012916.05 case is an SLD case (MVD = 150um). 
The top row, Figures 5(a) and 5(b), are the roughness and thickness for the 012816.05 case, and the bottom row, 
Figures 5(c) and 5(d), are the roughness and thickness for the 012916.05 case. 

Figure 5 demonstrates that while the icing event times are the same, the roughness and thickness topographies 
are very different. The magnitudes of the maximum roughness height and the maximum ice thicknesses are not 
noticeably different. However, because of the difference in stagnation point collection efficiency, caused by the 
larger droplet size, the ice roughness and thickness extends much further along the airfoil surface in the flow 
direction for the SLD case. 

Figures 6 and 7 present the spanwise average ice roughness and thickness measurements measured at each 
codebook vector of the SOM representation as a function of surface distance along the airfoils. Figure 6(a) presents 
the RMH variations for the Appendix C cases, while Figure 6(b) presents the RMH variations for the SLD cases. In 
addition to demonstrating the increasing roughness values as functions of time at any location, Figure 6 reinforces 
observations from Figure 5 regarding the ice thickness and roughness regions. That is, while the magnitudes of the 
maximum RMH values for any one case are similar, the roughness variations along the surface are very different 
between the Appendix C and SLD cases. Most importantly, the length of the glaze-ice plateau in the flow direction, 
as demonstrated in Figure 7, is significantly larger for the SLD cases than for the Appendix C cases. This extended 
plateau region suggest a longer liquid film region before film breakdown leading to the maximum roughness values 
of Figure 6 occurring further aft for the SLD cases than for the Appendix C cases. 

Figure 7(a) presents the average ice thickness values evaluated at each SOM codebook vector for the Appendix 
C cases, while Figure 7(b) presents the average ice thickness values for the SLD cases. Just as was the case for the 
roughness values, the maximum ice thickness values are not significantly different, but the distributions of ice along 
the surface of the Appendix C cases is much different than the distributions along the surface for the SLD cases. 
While the distributions are different, the most significant difference is that the nearly flat region near the attachment 
line (since the wing is swept) of the wing, which represents the glaze-ice plateau, is significantly wider for the SLD 
cases than for the Appendix C cases. 


B. Scaled Roughness and Thickness Temporal Variations 
Figures 8 and 9 present the airfoil maximum roughness (AMR) and attachment line thickness (No 47 ) 
evolutions in dimensionless form. Since the wings used in the study have a symmetric airfoil shape, the AMR 


values are defined as the average of the maximum RMH value on the left side of the wing in Figure 6 and the 
maximum RMH on the right side of the airfoil. 


+ RMB# 


max,Right max Left 


2 


RMH 


(19) 


In Figures 8 and 9, the AMR values and stagnation point thickness measurements are nondimensionalized by 27, 
while the ice accretion times, which are expressed as attachment-line fully-dense rime ice thicknesses, are 
nondimensionalized by 2ro as well. The dashed line shown in Figures 8 and 9 represent a C; value of 0.5 in Eqns. 
(10) and (12) as found from the unswept-wing roughness evolution investigations. Figures 8 and 9 demonstrate that 
the scaled AMR values as well as the scaled attachment-line thickness values for swept wings are increasing at rates 
nearly identical to the rates found on straight wings. 
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Figure 5. Surface Visualizations of the 100-second Cases: (a) IRT012816.05 Roughness Topography, (b) 
IRT012816.05 Thickness Topography, (c) IRT012916.05 Roughness Topography, (d) IRT012916.05 Thickness 
Topography (Top row: Appendix C, Bottom Row: SLD) 
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Figure 6. Surface Roughness Profiles along the NACA 0012 Surface Direction: (a) the Appendix C Cases and (b) 
the SLD Cases 
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Figure 7. Appendix C Thickness Profiles along the NACA 0012 Surface Direction: (a) the Appendix C Cases and 
(b) the SLD Cases 
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Investigating Figure 8 closely, for relative theoretical Rime ice thicknesses (Na,r/27ro) less than 0.1, the relative 
airfoil maximum roughness heights (AMR/2ro) are lower than the C; = 0.5 line. Figure 10 presents the airfoil 
maximum relative roughness values relative to the theoretical Rime ice thickness (AMR/ Na,z) versus the relative 
Rime ice thickness growth values. Figure 10 demonstrates that both the Appendix C and SLD cases exhibit an 
asymptotic behavior approaching a C; value around 0.6. The asymptotic behavior is reminiscent of a “first-order 
differential equation” effect, which would suggest an alteration to the temporal function of Eq. (10) to include a 
decaying exponential function as indicated in Eq. (20). 


3(r)=C, ( aes Wor (20) 


where to would be a time constant based on the cloud properties, flow characteristics, or possibly the liquid film 
dynamics. 

However, Figure 11 presents the relative maximum roughness values (AMR/ Na,r) versus A-cos(A), which is 
essentially the same abscissa as the relative, theoretical Rime ice thickness in Figure 10 without the influence of the 
cloud properties associated with the stagnation point collection efficiency, fo. Figure 11 demonstrates that the 
Appendix C and SLD asymptotic response collapse similarly to Figure 10. 

The fact that the data collapse similarly between Figure 10 with fo and Figure 11 without fo suggests that factors 
other than the stagnation collection efficiency of the cloud may be introducing a first-order transient response. That 
is, since the base cloud properties in the accumulation parameters are the same for all of the cases in this study 
resulting in an abscissa which is simply a scaled time function, the asymptotic behavior captured in Figures 10 and 
11 could be a product of wing internal conduction, conduction through the ice layer growing on the wing, or a real 
phenomenon associated with the roughness evolution physics. 

The transient artifact noted in Figures 10 and 11 may also be exhibited by the prior data for roughness evolution 
on straight wings. However, the ice accretion time progression measurements of McClain et al. [8] and McClain et 
al. [9] for straight 21-in. chord NACA 0012 airfoils have not been processed using the scaling factors employed in 
Figures 10 and 11. Understanding the true source of the first-order behavior, may require revising the straight-wing 
measurements; acquiring more measurements in the IRT for different cloud properties, specifically, different liquid 
water contents, freestream velocities; and acquiring roughness evolution measurements on airfoils and wings with a 
different construction from the solid aluminum wings employed in this study. 
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C. Surface Roughness Spatial Variation 

Figure 12 presents the roughness height variations of Figure 6 scaled by one-half the theoretical Rime ice 
thickness and plotted versus the local collection efficiency as defined in Eq. (16) for swept wing surfaces as was 
done similarly for Figure 2 for roughness on straight wings. In Figure 12, Figure 12(a) presents the Appendix C 
conditions, and Figure 12(b) presents the SLD conditions. Both Figures 12(a) and 12(b) indicate surface roughness 
distributions similar to those reported by McClain et al. [10] for straight wings. That is, both the Appendix C and 
SLD surface distributions exhibit regions near the glaze ice plateau with low values of roughness height (RMH/ Nar 
~ 0.1). Downstream of local surface collection efficiencies around 0.5, the RMH values rapidly increase to a 
maximum and then decay as the surface angles get closer and closer to 90° to the freestream. 

Also evident in Figure 12(a) and 12(b) is the first-order transient or asymptotic behavior noted in Figures (10) 
and (11). That is, at accretion event times less than 160 seconds, the measured roughness profiles are similar, but 
the magnitudes are increasing as the accretion time increases. After accretion times of 160 seconds, the measured 
roughness profiles exhibit essentially the same magnitudes. 

The gray dashed lines in Figure 12(a) and Figure 12(b) are the straight wing surface roughness variation trend 
presented in Eq. (12) using a liquid film breakdown collection efficiency (G,) of 0.45, which was the same film 
breakdown collection efficiency used for the unswept NACA 0012 cases of Figure 2, and a C; value of 0.5. Figure 
12(a) demonstrates that once the asymptotic conditions are approached, the swept wing Appendix C roughness 
distributions exhibit RMH values slightly above the straight wing curve. However, the surface distributions of the 
asymptotic Appendix C cases are captured by the straight wing distribution function. 

The asymptotic SLD results presented in Figure 12(b) do not collapse as well onto the unswept wing surface 
variation trend of Eq. (12) as well as the asymptotic Appendix C cases. As shown in Figure 12(b), the magnitudes 
of the RMH maximums is close to the unswept-wing maximum, but the SLD surface distributions appear wider, 
indicating a larger Weibull modulus or shape factor, and appear to be shifted downstream along the surface. This 
shift of the surface roughness profiles may be expected for cases with a much lower freezing fraction. While the 
SLD cases did have a slightly lower freezing fraction of 0.19 compared to 0.22 for the Appendix C cases, the 
difference was not expected to affect the surface roughness profiles as much as exhibited in Figure 12(b). 
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Figure 12. Airfoil Dimensionless Roughness Variations along Surface Related to the Local Collection Efficiency 


D. Final Observations Relative to Unswept Wing Roughness 

While the swept-wing SLD surface roughness distributions are slightly shifted from the unswept-wing surface 
roughness distribution and the Appendix C swept-wing roughness distributions, the swept wing ice thickness and 
roughness temporal evolution results of Figures 8 and 9 are almost identical to the temporal evolution of ice 
thickness and roughness on unswept wings. The asymptotic behaviors of the swept wing surface roughness 
distributions also exhibit the liquid-film dominated regions near the leading edge and the collection efficiency 
dominated regions near the impinging limits exhibited in the unswept wing results. Further, the asymptotic, swept- 
wing Appendix C cases exhibited nearly identical locations of liquid-film breakdown as was identified in the 
unswept-wing cases. While the sources of the first-order transient effects noted and the shift in the swept-wing SLD 
surface roughness profiles require further investigation, the results suggest that asymptotic swept wing and unswept 
wing ice roughness temporal and spatial variations may be modeled using: 


Weil B, —0.05, B p, —0.15,3) a 0 
+ 


RMH(N 4. 2B, Ba )=CNa.r Wei . tanh[10(8 », -a.}} (21) 


(0.875(2 ;, 0.15} By, —0.153) (2 2 
where 
Nar = 2%A-Bo cos(A) (22) 
B, = By cos(y, +ax)cos(A) (23) 
Bp = Bo cosly p +a)cos(A) (24) 


When the sweep angle, A, is 0°, Eqns. (21)-(24) reduce to the expressions provided in the Introduction section for 
unswept wings. 

An additional limitation on Eq. (21) relates to the effects of the stagnation point freezing fraction, n. The range 
of freezing fractions employed in this study and in the prior laser-scan ice roughness investigations is 0.19-0.25. 
Different freezing fractions are expected to affect the roughness surface distribution in three ways: 1) different 
freezing fractions may lead to different C; values, 2) different freezing fractions are expected to change the location 
of the liquid-film breakdown (fw), and 3) different freezing fractions may increase the width of the roughness 
distribution which could be accounted for by different values of Weibull modulus. Future efforts are required to 
determine the importance of and modifications required to account for the effects of different stagnation point 
freezing fractions on the surface roughness spatial distributions. 
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V. Conclusions 


As part of an extended study on the physics of ice roughness evolution and ice roughness surface distributions, 
the characteristics of ice roughness on swept wings were investigated using ice accretion time progressions for 
Appendix C and SLD conditions. Both Appendix C and SLD conditions were investigated using a 36-in. NACA 
0012 wings with 30°-sweep. The time and collection efficiency based scaling system presented by McClain et al. 
[10] for unswept wings was modified to consider wing sweep. Using the scaling modifications for sweep, the 
temporal and spatial roughness variations on the swept wing were compared to the previous roughness 
investigations for straight wings. The primary findings are: 


1) The temporal variations in maximum roughness height and attachment line ice thickness followed the 
straight wing results when the sweep-modified scaling parameters were employed. 


2) The temporal variations exhibited a first-order-differential-equation based asymptotic or capacitance-like 
behavior in time. While the effect may be a real phenomenon associated with the roughness formation 
process, airfoil internal conduction or ice layer internal conduction cannot be eliminated as contributing to 
the transient behavior observed. The prior unswept wing measurements from McClain et al. [8] and 
McClain et al. [9] must be revisited, and more tests may be required to determine the source of the 
asymptotic behavior in the scaled roughness results. 


3) Once the asymptotic time value was reached, the spatial variations exhibited on swept wings in Appendix C 
conditions were consistent with the spatial variations on straight wings. 


4) Once the asymptotic time value was reached, the spatial variations exhibited in SLD conditions were 
similar to the behavior found on straight wings, but the location of the glaze plateau was shifted further 
downstream and the roughness distribution along the airfoil surface was wider than exhibited on straight 
wings. The source of the spatial differences for the swept-wing SLD cases is not well understood. 


The results of the investigation suggest that the basic physics of roughness formation and evolution on swept 
wings are the same as those on straight wings and involve the competing phenomena associated with liquid-film 
region near the stagnation point or attachment line and the local collection efficiency downstream of the location of 
liquid-film breakdown. This observation is important because of the influence of the roughness elements on the 
formation of larger three-dimensional ice features formed on swept wings such as scallops. While the eventual large 
three-dimensional features on swept wings are very different from the primarily two-dimensional features on 
unswept wings following long ice accretion event times, both types of features begin as roughness that is evolving 
similarly and from the same mechanisms on straight and swept wings. 
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